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ABSTRACT 


21, INTRODUCTION 


The understanding of the instantaneous motions of turbulence is problematical owing to frequent unexpected 
changes. The central process in turbulent flows is convection induced by the instantaneous fluid velocity. 
Molecular diffusion fails to afford an impressive contribution to spatial transport at a high Reynolds number, and 
therefore convection controls the transference of momentum, chemical species, and enthalpy. Statistical measures 


hold well for the description of turbulent motion. 


A statistical description of turbulence comprises a probability distribution for stationary flows to determine the 
nature of turbulent velocity fluctuations. To study the statistical properties of turbulence Monin and Yaglom [1] 
enunciated the velocity structure functions (VSF) which are directly related to the probability density function 
(PDF) of the local dissipation rate. Therefore, an appropriate model for the PDF of the dissipation rate completely 


describes the statistics of turbulent velocity. 


The presence of intermittent behaviour of turbulent Lagrangian velocity statistics was discovered in various 
investigations on particle tracking [2, 3]. The massive progress in experimental procedures [4, 5] of measuring 
particle trajectories has caused a speedy development in the dynamics of tracer particles in incompressible flows. 
Some intriguing features of Lagrangian turbulence such as the notable role of coherent structures or almost 
singular structures compared to the Eulerian description have been revealed by numerical simulations [6], 


multifractal, and PDF-modeling [7,8]. 


In Arn*eodo et al. [9] an elegant collection of eight data sets from experiments and numerical simulations on 
turbulent velocity statistics along particle trajectories was carried out. The multifractal theory, which was extended 
to the dissipative scales and the Lagrangian domain, was seen to hold well for the studies in the intermittency of 
velocity statistics that were investigated. Intermittency attributes the multi-time correlations to the 


turbulence-rousing force. In effect, it was used in modelling the trajectory of a fluid particle by a multifractal 
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random walk [10] wherein long-time correlations and the incidence of large amplitude events at small scales 
dominate the motion. Based on this, the multifractal properties of the Lagrangian velocity may be explained. 
Lagrangian intermittency is associated with the nature, distribution, and time evolution of the dynamical structures 
entrenched in the flow. Such intermittencies were studied by using the velocity structure functions, which submit 
to the Kolmogorov similarity theory (KST). In fine, the central quantities of concern in the statistical theory of 


Lagrangian turbulence are the pdfs of the fluid particles [11, 12, 13, 14, 15]. 
“<2, LAGRANGIAN DEPICTION OF TURBULENCE 
The notion of temporal velocity increments along particle trajectories, 5)'(¢) =y'(t+¢)—v'(t) Which is the 


Lagrangian analogue of the Eulerian spatial increments of the velocity field can be introduced as an extension of 
the K41 phenomenology, originally advanced in the Eulerian framework. Here, v’ encodes the Lagrangian 
velocity, measured for each fluid parcel along its trajectory. The Statistics are assumed to be only dependent on the 
time increment, ¢ (i.e. stationary). The multi-scale temporal dynamics of Lagrangian temporal may be examined 
by observing the spanning statistical behaviours at different values of ¢. The multi-scale classification is largely 
realized by considering the dependence of the statistical moments of Lagrangian increments with the time 


increment, ¢ . These define the Lagrangian function structures: G+ (f)= ( dvi (oy ). 


2.1. Equations of motion 


Let X(t;v,y) be a transformation that maps the initial position X(t = 0; v, y) = y of a hypothetical tracer particle with 
initial velocity v onto its position X at later times ¢ > 0. The equations of motion of a single particle are of the form 


[12]: 


AX(GN,Y) _ yr. (1) 
dt ar UG; v.y) 
dU(t;v,y) : (2) 
= A : 
a (VY) 


where in (1) d/dt is the time derivative along the particle’s trajectory and the particle’s acceleration. In (2) d/dt is 


the Lagrangian or material time derivative that may relate to the partial derivatives of the Eulerian fields (see [16]) 


‘ 3 
By = (X(,1) +0 (KODA KO. . 
dt dt Ox, 
The acceleration, A in (2) can be expressed in terms of the Eulerian pressure and the velocity field estimated at the 
location of the particle. The particle’s position and velocity in the Eulerian sense, having the velocity field u(x, 1), 


related to the Lagrangian sense reads: 


dX 5 4 
—(t) =U, (1) =u, (X(),D) (4) 
dt 
and in a similar mode, the particle’s acceleration A(t) reads 
_dU,,. d’X,,.. (5) 
A(t= a =e (t) 
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In the Lagrangian mode, Navier-Stokes equations can be written in the form: 


; (6) 
dU, 1 Gp _., OU, 
at (1) pds (X@),0) + "an Be, (X(@),2), 
a 
ou, _9 (7) 
O*x, 
2.2. Lagrangian pdf. 


Let Xo = {Xo01, X02, X03} be Lagrangian coordinates and let fo be a reference time. Therefore x(t, Xo) denotes the 
position at time ¢ of the fluid particle that is at xo at time f [i.e. X' (fo, Xo) = Xo, where for each Eulerian variable U(x, 
t) say, the corresponding Lagrangian variable is denoted by the superscript (or subscript) L ]. From Pope [17] we 


get 
U'@,x,) =U[(@x,),¢] (8) 


Consequently, for the Eulerian velocity, x" is the solution to the initial value problem, 


Ox" (t,X,) _ 
a 


L 
x (t,X,) =X, 


U',x,), |" (9 a, b) 


Considering a constant (density) property flow [17, 18, 19], the Lagrangian joint pdf of the fluid particle properties 
at time ¢, conditional upon their properties at time fo, for the event 
{U'(,x,) =V, x"((4x,)) =x} (10) 


is given by 


g,(V.x;2|V,.x, )° (11) 


subject to the condition U'(to,Xo) = Vo. The relation (12) is analogous to the (one-point, one-time Eulerian) joint 
pdf of velocity g,(V; x, f). In general, the joint pdf of the events 
{U'(t.x, =V,, x'(t.x,) =x 37 =1,2....,N} (12) 


for N times, subject to the initial condition, Ut Xo) = Vo, is 


a (¥ 6Xtt, VAX. 


N-1?*"N-1? N-1?°° 


BV X30 |Viks) (13) 


“23, FLUID-PARTICLE VELOCITY 


The analysis of one component of the fluid-particle velocity U'(t) leads to a stochastic differential equation (SDE) 
that represents the Langevin equation, U(t). The crux of the velocity stochastic turbulence models is in the 


decomposition of a particle’s acceleration in the form 


“ = A(t) +7A(1)” (14) 
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where ¢A(t) and nA(t) encode the slow acceleration and fast acceleration respectively (see Sawford and Guest 
[20]. Whereas €A(t) evolves on the same timescale 1, as the velocity U" with 1, being of the order of k/e , nA(t) 


evolves on the timescale t, which is of the order of the Kolmogorov timescale 1,,. 


In a fully developed turbulence z, << 1, and therefore for the large scales, 


2 -(#) . (15) 
T, 


The generalized Langevin model for a velocity stochastic turbulence, in line with (14) is of the form [21, 22, 16] 


1 
q, -(a)4tu's.¢, (16) 
dt ns ’ dt 
In (16) above we let EA, =(a)+TU" and er i 
"dt 
wherein 
EA,(t)=(a,)(X(0,t) +7, (X(,1)| UO -(u,)(X(.1) |: (17) 


encodes the slow acceleration; 7; is the drift tensor which depends only on the location of the particle. The white 


noise that models the fast acceleration 7A,(t) is demonstrated by 


dW, 


7A (t) = D,(X(t),t) = (t) (18) 


where D is the diffusion tensor which also depends only on the location of the particle, and W(d) is a vector-valued 
Wiener process. For an insight into the diffusion tensor, let us consider a coupled stochastic process in velocity, u 
and acceleration a of the form (see Zamansky [23]): 


du, =adt, (19a,b) 
da, = Badt + D,dW, 


where dW; are the increments of the jth component of the Wiener process ( (dW;) =0; (dW.dW;) =dt6;); each 
of the vector B (the drift) and the diffusion tensor D depend on the vectors a and u. The tensor D can be 


decomposed into 


D =d6.+S.4+Q,  (dnotaconstant), (19c) 


where Sj; encodes a zero-trace symmetric tensor and Qj is an antisymmetric tensor. For statistical isotropy 
consideration of the acceleration, S;, must vanish but Q, need not be 0. The expressions for B; and Dj, are elegantly 
specified in [23]. Now, consider a stationary homogeneous isotropic turbulence, with zero mean velocity, 
turbulence intensity (TT) uw’, and T the Lagrangian integral time scale. A one-component Langevin equation, in U*, 
of the fluid-particle velocity U'(t) reads 


dU ‘(t) = UOT + Qu" IT AW) (19d) 
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with the rms fluid-particle velocity as u’. Note that U*(t) indicate the modelled particle properties for U'(t).The 


Lagrangian velocity auto-correlation RL(¢, s) given by 


(U'()U' (t+ 6)) (20) 
(UM) (UY t+O)))- 


R(t.o)= 


However, in stationary homogeneous turbulence, the Lagrangian velocity auto-correlation, R, (¢) depends only on 
the time lag,¢ , and is an even C’-function (see [24]). Thus, at the origin R'(¢=0)=0- The auto-correlation in 


(20) is the exponential function 


RQ) =e*" > R(C =0)=-1/T £0 > (21) 


and the time scale T satisfies 
P=[ Riojde. (22) 


We revert the equation (19) in which the Markov process U*(f) ‘seeks to model’ the temporal fluid particle 
velocity Uo. For a given initial condition at time tp, U*(fo) is a Gaussian random variable with zero mean and 
variance u’’; for t > to, U*(a) describes the Ornstein-Uhlenbeck (OU) process (the stationary random process; as 
the Wiener process is nowhere differentiable, the Langevin equation is, in a strict sense, only heuristic). As a 
constant drift term characterizes the Wiener process, the Ornstein—Uhlenbeck process holds well when it is 


dependent on the current value of the process. The latter is true in the present case. 


The intrinsic inadequacy is that while U'(t) is differentiable, U*(#) is not. Therefore the model fails a qualitative 
test if U*(f) is examined on an infinitesimal time scale (see Pope [17] and Minier [24]). However, a resulting 
brilliant analysis [17] using the behaviour of the autocorrelations function shows that much as an uninterestingly 


negative slope is recorded at the origin slope at very small times ¢/7 , owing to U*(¢) being not differentiable the 


exponential form provides a very reasonable approximation to the observed autocorrelations for larger times. 
<4, INTERMITTENCY 


Particle trajectories, under normal circumstances, submit to Gaussian statistics. Intermittency arises when 
deviations from Gaussian statistics get increasingly larger when considered at increasingly smaller scales of 
fluctuations. It is supposed that the small scales of turbulence are produced due to several interactions down the 
flow and are, at a reasonably large Reynolds number, independent of the large fluctuation scales and the turbulence 
production mechanism. The small scales depend uniquely on the mean rate of spectral energy transfer, which 
under stationary conditions amounts to the mean rate of dissipation of turbulence kinetic energy (s) , and the 


viscosity v. 


At suitably high Reynolds numbers Sawford and Yeung [30] indicated that the inertial sub-range, which is a 
portion of the equilibrium range, is detached from both the energy- containing and dissipation scales, and therefore 


independent of viscosity. It is, in effect, characterized solely by () . The dissipation sub-range is the portion of 
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the equilibrium range at the smallest scales where the viscosity is essential. Momentous as it is, velocity statistics 
at various temporal scales may depend on large scale forcing and boundary conditions, as evinced by the concept 
of universality. Admitting the central nature of the concept, Arn*eodo et al. [8] investigated intermittency and 
universality properties of velocity temporal fluctuations obtained from some laboratory data. The measuring of 
intermittency features of the velocity fluctuations is typically realized through the analysis of the velocity structure 


functions. 


Regarding moments of the velocity fluctuations, the Lagrangian Velocity Structure Functions (LVSF) of positive 


integer order p in the form 


(23) 


(jovi G)")= (e+ g)-vOl)=G""@? 


where i = x, y, z encode the velocity components along a single particle trajectory, the average is defined over the 
ensemble of trajectories, and the intermittency expresses itself in the anomalous scaling exponents, i, admitting 


the power law scaling [8, 25] 
(jov,(g)"): $*- (24) 
Giving to Kolmogorov (K41)[25] prediction 7*= p/3 with exponents that deviate considerably from the simple 


scaling prediction, especially for p > 3, where 1, <A“. The scaling exponents i, are determined by Taylor’s 
Pp 


expansion 
A, = Xa,(-1)"" p' 1k! (25) 


As indicated by Arn‘eodo et al. [8], the decay of the power law is imminent for scales ¢ _, ¢ , when dissipative 
effects dominate. The statistics of velocity fluctuations at changing time lag¢ can be quantitatively expressed by 


the logarithmic derivatives (see [26, 27, 28]) whose Taylor expansion is in the form 


(26) 


log ([0u,(6)|") = PAC)+ 2 A)+E AC) + 
where 4 (¢) encodes the pth-order cumulant of log| dv (|: For some derivatives Gi (¢) againstGn(¢) , we get 


dlogG"(6) (27) 


aA = 
(Ps) dlogG?(¢) 


Since 4 (p,¢) 18 ¢ -dependence, a scale-by-scale characterization of intermittency may be achieved. 


Considering the energy spectrum equivalence, the LVSF, and its higher-order extensions resident in the inertial 
sub-range have the Kolmogorov similarity applied to Lagrangian statistics defined along fluid particle trajectories 


through the flow in the form 
(via+ g) -vi@’) : (ey Cece T) (28) 
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and 
(pergy-vi@ol): wi, = 7,) (29) 


The dissipation (Kolmogorov) scales are such that length @ =(v° He)", velocity Vv, = (v(e))"" and time 


Ce (v(e))" > where Uencodes the viscosity term. 


The deviations from Kolmogorov similarity theory are considered a culprit in the event of intermittency. At small 
scales the dissipation rate of turbulence kinetic energy is not smooth, being an extremely variable quantity. In 
consequence, the cascade process is not sufficiently described by (¢) alone, and higher moments of the dissipation 


rate are important (see Sawford and Yeung [29]). 


In temporal (Lagrangian) terms a specific value of the local dissipation rate is 
_ 1 +o 7 ’ (30) 
E.(t)= ct é'(t')dt' » 


where «'(t) = e(x'(t), f) is the rate of dissipation along a fluid particle trajectory. The point-wise values of the 


dissipation rate are obtained as € — 0. 


The Refined Similarity Hypothesis (RHS) admits writing equations (28) and (29) in the form 


(ve+o)-vi@ 7 le,): ae - (aca 7y (31) 


and 


(wE+g—v®) 


: |e.) [eu 6 aa (32) 

The dissipation sub-range RSH scaling was tested using (32) for (¢ = i) (see [29, 30] for details) by using the 
conditional structure functions. It was seen that for large values of é./(e) the conditional structure functions 
approach the dissipation sub-range RSH scaling but there is a strong departure from RSH scaling as E. / (e) — 0. 
As ¢ — 0 the corresponding non-dimensional conditional moments of the acceleration are recovered. 

== 5. CONCLUSION 


The most classical problems of Lagrangian turbulence relate to the issue of particle dissipation, with two salient 
aspects: (a) diffusion of single particles from a point source (typically, Taylor problem) and (b) the relative 
dispersion of pairs of particles (the Richardson problem). The Lagrangian velocity may be explained by modelling 
the trajectory of a fluid particle in a multifractal random walk in which long-time correlations and the incidence of 
large amplitude events at small scales dominate the motion. The associated nature, distribution, and time evolution 
of the dynamical structures entrenched in the particle flow introduce intermittencies in the flow regime. The LVSF, 
which submits to the Kolmogorov similarity theory (KST) was applied in studying such intermittencies. Some 
findings suggest that deviations from the Kolmogorov similarity theory are considered a culprit in the event of 


intermittency. 
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